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Abstract 

A QCD sum rule calculation of Balistky and Ji on the spin content of 
the nucleon is done with a different approach to the evaluation of the bilocal 
contributions and to the extraction of the nucleon pole residues. The result 
obtained is much more numerically stable which puts their conclusion that 
more than half of the nucleon spin is carried by gluons on firmer ground. 

The question of the spin distribution of the nucleon among its constituents 
has recently received much attention and it is now realized that a substantial 
fraction of the nucleon spin is carried by the quark orbital and gluon angular 
momenta [|TJ 

~ = \^ (/*') + L q (f) + Jq {^) (1) 

Where fj 2 is the scale at which the relevant operators are renormalized (we 
shall take here /i 2 ~ lGeV 2 and omit it from the notation). The first term in 
eq.(|l|) has been measured in polarized deep inelastic scattering 0. The second 
and third terms represent quark orbital and gluon contributions respectively. 

Ji H has recently shown that the total quark (and hence gluon) contribution 
to the nucleon spin is measurable through virtual Compton scattering in a special 
kinematic region where single quark scattering dominates via the study of off 
forward parton distributions (OFPD) whose second moment yields the quark 
contribution. 

It was further shown by Vanderhaegen, Guichon and Guidal that hard 
electroproduction of photons and mesons on the nucleon also gives access to the 
OFPD. 

These authors also suggest that preliminary measurements can soon start at 
HERMES and CEBAF which can shed light on the spin content of the proton. 
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An explicitly gauge invariant decomposition of the angular momentum oper- 
ator in QCD was first given by Ji @ 
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This has subsequently been used by Balitsky and Ji to study the three point 
function 



Wf a (p) = / exp(ipx) exp(iqz) (0 | T rj(x) fj(0) M^ a (z) | 0) dx dz (3) 

y limg— >0 ' 



where rj is the nucleon field 



rj(x) = e y >"c7 a v?) 7s la <f 

M M7a = T m x v — T^ u x a is the angular momentum density 
with 



T af) = T aP + T a/3 = ^ (a ■ ^ 0)^ + _ p<*H pP ^ 

The energy momentum tensor of QCD written as the sum of its quark and 
duon parts ((a/3) means symmetrization in the indices) 



W^ a = W(p 2 )(2ip' M Y^ a ) + other tensor structures (4) 

W(p 2 ) contain nucleon double and single pole contribution as well as a non sin- 
gular contribution of the continuum 

W(p 2 ) = / gA % + 2 d +•■■ (5) 
(p 2 - m^) 2 (p 2 - mfv) 

A at is the coupling of the nucleon to the current rj 



(Oh | j\T(p)> = Ajy w(p) (6) 
Using a Ward identity the 3-point function(|]) is rewritten as 



Wj?" = dx dz z v exp(ipx) exp(iqy) (0 | T rj(x) fj(0) O a (z) | 0) - v «-> a 

limg^O J 

O a (z) = i; f F a ? 7 ^(z) (7) 
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In the deep Euclidean region W can be calculated in QCD 
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The terms in In 



signal the breakdown of the validity of the OPE for 



q — > 0, they have to be replaced by bilocal contributions ||, i.e. one first 
expands T rj(x)r)(0) = c n( x )O n for small values of x, the operators O n then 
combine with O a (z) in a series of two point functions the leading terms in the 

asymptotic behaviour of which reproduce the terms in In ^— ^\ and which must 

be extrapolated to q 2 = 0. Expression (||) must then be replaced by 
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and the relevant operators 
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In the deep Euclidean region 
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The method followed here to extract Ho{q 2 — 0) differs from the one used 
in ref||. n (s) is an analytic function in the complex s-plane except for cut 
along the positive real axis. Consider the integral over the closed contour in the 
complex s-plane consisting of a large circle of radius Ri and two straight lines 
immediately above and below the cut which run from threshold to R\. By virtue 
of Cauchy's theorem 




= m /2 n (0) (13) 

Where m! is a mass parameter. The integral over the real axis corresponds to 
hadronic intermediate states with quantum numbers J pc = 1~ + , as m' 2 is varied 
from threshold to R this integral changes sign which implies that it vanishes for 
some intermediate value of m' 2 which is expected to be close to the value of s for 
which Imllo(s) reaches its maximum. Experimentally the E857 collaboration 
[U has recently observed a broad bump in the invariant mass squared of the r]n 

system produced in the reaction n~p > r]7i~p around 2.0GeV 2 . So if we take 

m' 2 ~ 2.0GeV 2 we expect the integral over the real axis in eq. (|i~3|) to be negligible. 
In the integral over the circle Ho is well approximated by IIq except possibly 
for a small region near the real axis. Eqs. fll~2]) and (|3~3|) then yield 

2 a s R\ ( Ri 1\ 1 ,a s F 2 f R x R x \ 8 2 1 

n ° (P = 0) " 60^ - 2 ) + T2 ( ^T } - L VJ ~ ^ as{{qq) ^ 

(14) 

In a similar fashion the asymptotic behaviour of IIi(s) 

n? CD (*) = -| ^(qq 2 ) In + (gqFaqqq) + ■ ■ ■ (15) 

yields 

ntf = o) = H^( ? V> (^-"4) +\4m«) 2 ) (w) 

54 7r \m u fi J 6m' A 

where we have parameterized 

{gqFaqqq) = -ml((qq) 2 } (17) 
In the equations above Ri ~ 4GeV 2 and ml ~ .65GeV 2 are reasonable estimates 

E 
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We are now able to extract the residue J g . This is done by considering the 
Laplace type integral J dt exp (— -g^) W(t = p 2 ) in the complex t-plane 
over a closed contour consisting of a circle of radius R and two straight lines 
above and below the cut which run from threshold to R. M 2 is the usual 'Borel 
mass' parameter. The exponential provides convenient damping of the integral 
over the continuum which we expect to be small for an appropriate choice of M 2 . 
If this is done we get a relation between the residues of W and the integral of 
\yQC~D on c j rc } e Q f radius R 
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X N is obtained from a study of the two point function 



/ dx exp(iqx) (0 | Tr/(x) r/(0) | 0) 



(27r) 4 A 2 ,expf-^ = -M 6 dxx 2 exp (-x) - — (^)M 2 rfx exp (-x) 
V M 2 y 4 Jo 2 7T J 

QO 

+y^ 4 (fe) 2 ) (19) 

The choice of M 2 in the equation above and the coherence of the method 
is dictated by stability considerations, if there are values of M 2 small enough 
to provide adequate damping of the continuum and large enough to justify the 
neglect of the contributions of higher order condensates in the OPE, this should 
show up in the linear behaviour of the r.h.s. of eq.(|l8|) in some intermediate 
range of M 2 . The standard value § (^) = 0.012GeV 2 is used, for the value of 
the condensate ((qq) 2 ) the choice (qq) 2 (vacuum saturation hypothesis) is usually 
made but as it seems too stringent an assumption we take ({qq) 2 ) = (3(qq) 2 
the dependance on (3 is weak because it appears both in the numerator and in 
the denominator of the r.h.s. of eq . (|i~8l) . In fig.l the r.h.s. of eq . (|T8|) is plotted as 
a function of M 2 for (3=1. It is seen that the dependance on M 2 is essentially 
linear in the interval .7GeV 2 < M 2 < 1.2GeV 2 with a small value of the slope 
{d ~ 0). 
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The result is 



J g ~ 0.35 for p = 1 , J g ~ 0.30 for (3 = 3 (20) 

Balitsky and Ji || had reached the same conclusion, they had however chosen 
to eliminate the contribution of the nucleon single pole by multiplying W by 
(t — m 2 N ) before performing the Borel transformation and used vector meson 
dominance to estimate no,i(0). This has the unfortunate effect of destabilizing 
the calculation as it can be seen that the variation of their value for J g with M 2 
is very rapid. 

This can be seen explicitely from the expression obtained following the method 
of ref.[||]. Instead of eq.(pT8|) one gets 




stability here requires the r.h.s. of eq.([H]) to be constant in an intermediate range 
of M 2 . As can be seen from fig. 2 no such thing happens. This confirms what has 
been stated above: elimination of the simple pole contribution by multipication by 
(t — rrvfy) before perfoming the Borel mult if or mat ion destabilizes the calculation. 
The main difference between the approach advocated here and the one of ref. ||, 
in addition to the use of eq . (|13"D instead of vector meson dominance to evaluate 
the bilocal contributions no,i(0), consists in not multiplying by (t — m 2 N ) in order 
to restore numerical stability. 
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Figure Captions: 



fig. 1 : The r.h.s.of eq . (^8|) plotted against M 2 . The linear variation of 
J g + c'M 2 is apparent in the range .7GeV 2 < M 2 < 1.2GeV 2 

fig. 2: The r.h.s. of eq.(|21|) plotted against M 2 . J g which ought to be constant 
in some range of M 2 varies rapidly. 
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